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Abstract 

Hard single diffractive processes are studied within the framework of the 
triple-Pomeron approximation. Using a Pomeron structure function motivated 
by Regge-theory we obtain parton distribution functions which do not obey mo- 
mentum sum rule. Based on Regge- factorization cross sections for hard diffrac- 
tion are calculated. Furthermore, the model is applied to hard diffractive particle 
production in photoproduction and in pp interactions. 



1 Introduction 



A significant fraction of events in high energy hadron-hadron interactions and 
in photoproduction is characterized by a rapidity gap between a quasi-elastically 
scattered primary hadron and a multiparticle final state. These single diffractive 
interactions can be understood in terms of Pomeron exchange. 

As suggested by Ingelman and Schlein [|IJ there is now - in addition to the 
well-investigated soft diffractive particle production - experimental evidence for 
single diffractive jet production both in pp and ep interactions 0, [| f|, [^, [| [7|]. 
Characteristic features of these jets are very similar to the ones produced in 
nondiffr active events. This observation suggests that it might be possible to 
apply models based on perturbative QCD to describe hard diffractive scattering 
processes between hadrons or photons and Pomerons. 

It is useful to distinguish between noncoherent diffraction and processes where 
the Pomeron enters the hard scattering as a whole. Models based on noncoherent 
interactions usually assume the Pomeron to be a partonic object, the parton- 
momenta obeying parton distribution functions (PDFs) Jj], §, |[ |l(], [□], [12| . Ac- 
cording to the assumptions made, these PDFs strongly differ in the partons con- 
sidered, in their shapes, and in their normalizations. In addition one can consider 
interactions where almost all the Pomeron momentum enters the hard scattering. 
First signs for a "super hard" Pomeron structure in pp interactions were reported 
by the UA8-Collaboration 0]. These effects may manifest itself in a breakdown 
of factorization and would lead to a delta-function like term in the Pomeron 
structure function. They can be explained by mechanisms called "coherent hard 
diffraction" |13|] or by means of a direct Pomeron-quark coupling [^, [14|, |15 



However the experimental information on diffractive interactions containing jets 
in the final state is still limited and does not allow to draw definite conclusions 
on the Pomeron structure. 

In the present paper we discuss hard diffraction in the framework of the two- 



component Dual Parton Model [16, 17, 18, 19]. We use an ansatz developed by 



Capella et al. |2T| which is based on Regge-theory to obtain the Pomeron 
PDFs. Further assumptions on the normalization of the quark distributions are 
not necessary since this is given by Ff- Features of particle production in hard 
diffractive photoproduction and in hard diffractive pp interactions are investigated 
and compared to data [[| [5J. 

The paper is organized as follows. In Sect. 2 we describe the way hard diffrac- 
tion is treated in our model. We especially focus on the PDFs of the Pomeron. 
Using Monte Carlo realizations of the model, we discuss hard diffractive particle 
production and calculate hard diffractive cross sections in Sect. 3. A summary is 
given in Sect. 4. 
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2 Description of the model 

2.1 The triple— Pomeron approximation 

Within Regge-theory, high-mass single diffractive processes are understood by 
means of triple-Pomeron exchange (see Fig. p]). Using Reggeon field theory |p2[ 
this graph is treated as an usual Feynman graph. The Pomerons shown in Fig. [I] 
are defined by propagators £ip(£)(s/s ) aip ^' ) an d cannot be considered as parti- 
cles In lowest order the differential diffractive cross section is given by the 
unitarity cut (Fig. P 

(1) 

With M we denote the mass of the diffractively excited system. g^,g^2, and 
r 3IP are the various couplings as shown in Fig. [I]. is the usual signature 

factor 

sm.(irajp{t)) 

ajp(t) = 1 + A + a'jp(0)t is the Pomeron trajectory with the intercept 1 + A, and 



Sq = lGeV 2 . Introducing an effective Pomeron-particle cross section [23 



T 3JP (t 0) ( M 2 \ aip{0M 
< o T(M 2 ,t)=^(0)LJ^l ^_ , (3 ) 

Eq. (H) can be interpreted as a product of a Pomeron flux factor and this cross 
section 

To estimate the contribution of hard diffraction (an example of such a process is 
shown in Fig. |2|) to the single diffractive cross section, can be replaced by the 
hard Pomeron-hadron/photon cross section af . We obtain af applying lowest 
order perturbative QCD, i.e. in case of hadron or resolved photon interactions 

of = E t^t t dx * I di ^^fa(^Q 2 )f^Q 2 )Q( P ^- P r off ) 

(5) 

and for direct photon-Pomeron interactions 

<ir = £ / W di ^^fU^ Q 2 ) &(p± - vT fS )- (6) 
j,k,i Jo J dt 
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We sum over all parton configurations and integrate having a lower cut-off in 
transverse momentum p c ^ to ^ . Here fl(x, Q 2 ) are the PDFs of the hadron/photon 
whereas with fjp{x,Q 2 ) we introduce distribution functions of partons inside 
the Pomeron. A detailed discussion of the Pomeron-PDFs will be given in the 
following section. 

The free parameters occurring in Eq. ([]]), such as the proton-Pomeron cou- 
pling constant and the intercept of the Pomeron-trajectory, are obtained within 
the two-component Dual Parton Model by fits to data on total, elastic, and 
diffractive cross sections. We refer to fl8], |19L p4l Ef| for further details. 



2.2 Parton distributions in the Pomeron 



As shown in |9|, 14), [12], |2l| one can derive a Pomeron structure function 



ip 



by relating it to the total cross section of virtual photon-proton diffractive deep 
inelastic scattering. Similar to the way the proton structure function F% can be 



related to aU one gets 



X 



Q 1 



Q 2 + M 2 



(7) 



Convoluting Ff with a Pomeron-flux factor, one obtains the structure function 
of diffractive dissociation. Using Regge factorization it is possible to calculate 
Ff from the deuteron structure function F$. The complete formalism is given in 



Ref. [2(J, [Hp. Here we only want to give the main formulas which are necessary 
to understand the PDFs of the Pomeron we will use afterwards. The proton 



structure function is parametrized at moderate values of Q (Q < 5 GeV ) [ 20 1 



F%(x, Q 2 ) = S(x, Q 2 ) + V(x, Q 2 



Ax~^ Q2 Hl 



x 



,n(Q 2 )+4 



Q- 



1+A(Q 2 ) 



-B x l - aR (l-x) n[Q2) 



^Q 2 + aJ 

Q 2 V R 



Q 2 + b, 



(8) 



with 



A(Q 2 



A 1 + 



2Q' 



Q 2 + d / 



n{Q 2 



3 
2 



1 + 



Q 2 



Q 2 + c / 



(9) 



We again refer to |20|, |2]J for the exact values of the parameters entering the 
expressions and a detailed discussion of the Q 2 -dependent intercept. The Pome- 
ron structure function can be related to the deuteron structure function by the 
following substitutions 



F*(x,Q 2 ,t) 



F*(x,Q 2 ;A 



e(t)A, B - f(t)B,n(Q 2 ) - n(Q 2 



(10) 
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e and / are ratios of coupling constants. The t-dependence of the Pomeron- 



PDF is completely given by the ratios e and /. We use e(0) = 3/(0) = 0.1 |f2T 
and an exponential dependence on t with the slope b = 0.5 GeV~ 2 ||23|. [19|, pp . 
The substitution n(Q 2 ) — > n(Q 2 ) — 2, for example, is due to the similarity of the 
valence quark distribution in the Pomeron and a meson. The first term in Eq. (|8|) 
determines the sea-quark distribution whereas the last term is responsible for the 



behaviour of the valence quark distribution. Using Eq. fllOD , and the definition 
F?(x, Q\ f) = £ 4*(/t <5 2 ; t) + &(x, Q 2 ; t)) (11) 
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where e q are the corresponding quark-charges, we obtain at moderate Q 2 -values 
for the Pomeron PDFs 

x fip( x iQ ]t) = x fw( x iQ \ t) = x fip( x i Q it) = x fw( x iQ it) 
= ^ e (t)S(x,Q 2 ) + ^f(t)V(x,Q 2 ) 

xf^(x,Q 2 ;t) = xf^(x,Q 2 ;t) = ^e(t)S(x,Q 2 ). (12) 

For simplicity, we assume a SU(3)-symmetrical sea of light quarks. The charm- 
quark distribution of the Pomeron is suppressed at the Q 2 -values considered so 
far. In contrast to the nucleon it is not possible to determine the normalization of 
the gluon distribution in the Pomeron by momentum sum rule. We are therefore 
free to choose a specific form and normalization. With respect to the experimental 
observations which favor a relatively "hard" structure we use 

x&(x,Ql;t) = K(C&,t)x(l-x) (13) 

where K(Q^, t) depends on the normalization chosen. In the following we take the 
normalization of the gluon distribution according to the scaling factor resulting 
for the quark distribution from Eq. (|T0"D. 

As already mentioned all distributions obtained are limited to Q 2 < 5 GeV 2 . 
In order to get PDFs of the Pomeron at higher Q 2 -values we take ( |I"2"D and (p~3|) 
as input distributions for a QCD evolution in the leading logarithmic approxima- 
tion |T2|. The evolution was done using the code of |27|, with Ql = 2 GeV 2 . 



The result is shown in Fig. ^]and[|, where we have plotted xfjp,xf^>, and xfjp for 
different values of Q 2 . Within this formalism, the normalization of the Pomeron 
PDFs is given by 

N q = Y, tdxxfl(x,Q 2 ,t), N g = f'dxxfU^Q^t) (14) 

which is shown for t — GeV 2 in Fig. [5[ The hardness of the u- and (i-quark 
distributions at low values of Q 2 is governed by the valence part of (0), i.e. 
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for x > 0.5 we are mainly dealing with a \/x{l — x) n( ^ 2 )~ 2 -behaviour, n(Q 2 ) 
being about 2.0. ..3.0 (c = 3.55 p0||). We note that this behaviour is similar 
to the Q ^independent predictions of ||, 0, whereas at high Q 2 -values our 
quark-distributions become softened due to the sea-quark contribution and the 
Q CD-evolution. 

Furthermore, we show in Fig. |6| the structure function for different de- 
scales. Again, the flat shape of is determined by the valence quark distribu- 
tions. 



Particle production in hard diffractive inter- 
actions 



3.1 Sampling of hard single diffractive events 

The Monte-Carlo treatment of soft single diffractive hadron-hadron interactions 
within the two-component DPM is described in ||T7| , 18fl . The generation of single 



diffractive events in photoproduction will be discussed elsewhere Here we 

want to focus on the way the existing models have been extended to diffractive 
jet-production. This has been done for pp collisions using an extension to the 



DTUJET-93 code [30] and for pp as well as 7p collisions using the DTUJET- 



PHOJET code §1 



The Monte-Carlo implementation is similar to that of usual hard scattering 
processes (i.e. processes involving high p t ) between hadrons or between hadrons 
and photons. The main differences are: (i) the interaction is boosted to the 
rest system of the diffractively excited "blob" , the CM energy is therefore given 
by the diffractive mass, and (ii) one hadron is replaced by a Pomeron with a 
virtuality t. The momenta of the partons entering the hard 2 — > 2 scattering 
process are obtained using conventional PDFs of hadrons and photons 33[ [33 



and the PDFs of the Pomeron introduced in the previous section. Initial state 
radiation which significantly modifies the multiparticle final state of interactions 
involving high transverse momenta has been implemented. The chain system to 
be hadronized using JETSET 7.3 |35|, |36| is determined by the color flow taking 
cross sections for different color flow diagrams into consideration |J7|, 3S|. 

Parts of the MC-realization of hard diffraction are technically similar to the 
ones described in [3U, iO"], which were so far mainly used to understand the under- 
lying interactions. However, we would like to emphasize that our starting points 
are completely different. Our investigation is based on the two-component DPM 
which treats soft and hard scattering processes in an unified manner. The free 
parameters are fixed by fits to cross section data. Similar to [^, f|(| we assume 
hard diffraction to be based on a partonic structure of the Pomeron but we start 
from a Pomeron structure function which is completely determined by fits to 
data on the proton structure function and a ratio of coupling constants which 



5 



follows from the model. There is no further freedom in choosing a specific x— and 
Q 2 -dependence of the quark distribution inside the Pomeron. The normaliza- 
tions of the Pomeron PDFs (and therefore the hard diffractive cross sections) are 
obtained from the model rather than imposing additional assumptions |39|, |40| . 



3.2 Hard diffraction in photoproduction 

Recent measurements at the electron proton collider HERA at DESY have shown 
that a substantial part of minimum bias photoproduction || |7j and deep inelastic 
scattering j|, || events exhibits diffractive features similar to hadron-hadron scat- 
tering. First distributions of so called rapidity gap events have been published 
[El F7], but the data are not yet corrected for acceptance and do not allow absolute 
comparisons. We will apply the model developed to investigate diffraction dis- 
sociation of photons in quasi-real photon-hadron scattering and we understand 
that more data will be published soon. 

The simulation of diffractive events in photoproduction can be done similar to 
diffraction in hadron-hadron scattering substituting the hadron-Pomeron scat- 
tering subprocess by photon-Pomeron scattering. In the calculation of the total 
diffractive cross section, the absorptive corrections due to multiple photon-hadron 
scattering are taken into account |19| . For the simulation of photon-Pomeron 
scattering we neglect unitarity corrections (e.g. multiple photon-Pomeron scat- 
tering) which become important only for very high diffractive masses^. Thus, the 
simulation of hard photon-Pomeron scattering follows directly from Eq. (||]||). 
The Pomeron flux is calculated using Eq. ([]]). 

To calculate ep photoproduction cross sections, the flux of quasi-real photons 
has to be estimated. Here we apply the improved Weizsacker-Williams approxi- 



mation 41 



A*M = f (i+MmSp-!^ (15) 



■ 1 1 tin 



with the kinematical cuts 

0.25 < y < 0.7 Q 2 mm = 3 • lO^GeV 2 Q 2 max = lO^GeV 2 . (16) 

y and Q 2 denote the energy fraction taken by the photon from the electron and 
the photon virtuality, respectively. 

Using our Monte Carlo program, complete hadronic final states have been 
generated and analyzed. According to the experimental conditions only events 
passing the r] max -cut |42|, § have been accepted for the further analysis. In ad- 
dition, a few less restrictive cuts have been applied to match the experimental 
selection procedure described in ||. 

In Fig. |7| we show the transverse momentum distribution of charged particles 
for the selected rapidity gap events in the pseudorapidity range —1.5 < r\i ah < 1.5. 

* High diffractive masses are suppressed by applying the experimental cuts. 
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The absolute cross sections obtained with the model are given. Together with 
the predictions we show uncorrected data of the HI Collaboration || on charged 
tracks. It is expected that the shape of this distribution can be compared with 
calculations for transverse momenta higher than 1 GeV/c |43[]. Note that the data 
are scaled in order to compare their shape to our calculation. The systematic 
difference between the model and the calculation at low p± can be qualitatively 
explained by the p^-dependent experimental acceptance in this region ||43|| . In 
addition, the contribution from direct photon processes is shown. 

Using a cone jet algorithm similar to the one used in || the transverse energy 
and pseudorapidity distributions of jets with E t > 4 GeV are calculated. In Fig. || 
we show the transverse energy distribution of jets together with uncorrected HI 
data 0. Since the so called jet-pedestal is small in diffractive events |5|, |7|, (43 



the transverse jet distribution should not be drastically influenced by acceptance 
effects. In Fig. |S] the corresponding pseudorapidity distribution of jets is compared 
to HI data. In both figures, the HI data are again scaled to compare their shape 
with the model predictions. 



3.3 Hard diffractive proton— antiproton interactions 



3.3.1 Cross sections for hard single diffraction 

An estimation of cross sections for hard diffractive events containing jets in the 
final state can be obtained using Regge-factorization (Eq. (f|)) together with the 
hard Pomeron-particle cross sections given in Eq. (|5|) and (0). However the cross 
sections depend strongly on the partonic cut-off in transverse momentum which 
enters the integrations. In addition, further uncertainties arise from choosing a 
certain scale in ([5]) and (||) and from the definition of a diffractive event itself, 
i.e. from the t- and M 2 -ranges the differential hard diffractive cross section has 
to be integrated over. 

Reliable predictions for cross sections in hard diffraction can therefore only 
be given for a certain experimental set-up taking into account all kinematical 
cuts applied and jet rates based on jet-finding algorithms which were used to 
obtain the experimental results. This has been done for the experimental set-up 
of UA8 || with a jet-finding algorithm which will be described further below. 
We calculate the ratio of the hard diffractive cross section and the total diffrac- 
tive cross section using the MRS DO' set for the parton distributions in 
the proton/antiproton and g p jp(0) = 6.2 Vmb, T 3JP (t = 0) = 0.08 \Jmb GeV 2 , 
0.25 GeV" 2 and Q 2 



A = 0.078, a'(0) 

spond to the UA8-cuts we multiply these values by the fraction of those events 



BJI. To obtain ratios which corre- 



which contain at least two jets of a transverse energy E{ > 
pseudorapidity \v\ < 2 in their final state. Experimental data for these ratios in 



preliminary form were given in [45] , final values will be available soon |^6| and we 



understand, that our calculations shown in Fig. 10 will be consistent with these 
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data |47| . 

3.3.2 Hard diffractive particle production 

Hard diffractive proton-antiproton interactions were recently investigated by 
UA8-Collaboration at a CM energy of 630 GeV. Only limited information 
on absolute hard diffractive cross sections is available so far j5| (see previous 
section). It is therefore not possible to compare our model directly to these 
data. However, since our model is able to describe data on hard diffraction in 
photoproduction rather well it might be worthwhile to show also predictions for 
hard diffractive pp-interactions. In particular we apply the same cuts to the fi- 
nal state as they were used in the experiments ||. Jets are identified using a 
cone-algorithm in the r] — 0-plane, with being the azimuthal angle and |?7| < 2. 
If a jet with E t > 8 GeV is found the search for its axis is iterated calculating 
i^-weighted sums over cells within an unit-cone radius. According to the UA8- 
data all distributions discussed in this section are normalized to unit area. Our 
results are obtained with the PDF set MRS DO' for the parton distributions in 
the proton/antiproton. 

As mentioned in [[| a variable sensitive to the partonic structure of the Po- 
meron could be Xj e t. It is defined as the longitudinal momentum component of 
a jet normalized to its maximum value in the Pomeron-antiproton CM system. 



In Fig. 11 we present predictions of our model for momentum fractions of the 
quasi-elastically scattered proton x p between 0.92 and 0.94. As stated by UA8 || 
the data were obtained with essentially full acceptance at positive x^-values, i.e. 
it can be expected that our calculations agree well with the data also after the 
corrections have been applied. 



In Fig. [12] we show the distribution of the jet-pseudorapidity in the antiproton- 
proton CM system again for the same x p -bin. The tail at positive //-values is 
influenced mainly by the structure of the Pomeron. 

A variable which may indicate whether there is a "super hard" Pomeron 
structure is X2j e ts - the longitudinal momentum of a two-jet system, again nor- 
malized to its maximum value in the Pomeron-antiproton system. The results of 
our calculation are plotted in Fig. ^ showing the contribution from quarks and 
gluons of the Pomeron-PDF separately. Due to the shape of the valence quark 
distribution (see Sect. j2.2| ) the quarks mainly contribute to higher values of xij ets 
whereas the gluon contribution is peaked around 0.2. Using our Pomeron-PDF 
we therefore obtain a significant fraction of events containing 2-jet systems with 
%2jets > 0.7. Considering the reasonable agreement found with the jp data the 
question to which extent a "super hard" contribution is still necessary within our 
model cannot be answered before corrected measurements are available. 
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4 Summary and conclusions 



Hard diffraction in hadron-hadron collisions and in photoproduction has been 
investigated in the framework of the two-component DPM. Since soft diffrac- 
tive interactions at collider energies are well-described in terms of Pomeron ex- 
change processes, diffractive jet-production may provide new information on the 
Pomeron-structure 

In the present paper we study features of hard diffractive particle produc- 
tion treating the Pomeron as a partonic object. PDFs of the PomeronjJ] are ob- 
tained from a Pomeron structure function motivated by Regge-theory. Whereas 
the quark distributions follow directly from the parametrization of the Pomeron 
structure function, it is possible to choose an ansatz for the gluon distribution 
which accounts for the experimentally observed hard Pomeron structure. The 
normalizations of the quark distributions of the Pomeron are determined by the 
scaling factor relating the Pomeron structure function to the deuteron structure 
function. The Pomeron PDFs are evolved to high values of Q 2 applying leading 
logarithmic QCD evolution equations. Using the triple-Pomeron approximation, 
we are able to give predictions on absolute cross sections and distributions. 

We demonstrate that our model and the Pomeron PDFs used are in reason- 
able agreement with presently available data on hard diffraction in 7^ collisions 
from the HERA-collider. Since there are no absolute distributions published the 
comparison was restricted to the shape. 

Cross sections for hard diffraction depend strongly on kinematical cuts and 
on assumptions defining an event as being produced diffract ively. Predictions can 
therefore only be given taking a specific experimental set-up into account. 

A comparison of hard diffractive particle production in hadron-hadron inter- 
actions to data is presently limited to data published by the UA8-Collaboration 
which do not allow absolute comparisons. Nevertheless, the results obtained 
within our model look promising and may explain the main features of these pro- 
cesses. We are not yet able to draw any conclusions concerning a possible "super 
hard" Pomeron structure. Further investigations will be necessary as soon as 
more data become available. 
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Figure Captions 



1. A cut triple-Pomeron graph (a) describes multiparticle final states charac- 
terized by a rapidity gap (b). 

2. Example for a hard diffractive scattering process. A gluon of the Pomeron 
undergoes a hard scattering with a gluon of the lower particle labeled with 
2. 

3. u- and s-quark distributions in the Pomeron are plotted for different values 
of Q 2 . The initial distribution of the QCD-evolution is shown for Ql = 
2 GeV 2 . 

4. Gluon distribution in the Pomeron shown for different values of Q 2 . 

5. Normalization of the quark and gluon distribution function of the Pomeron. 

6. Pomeron structure function F^(x, Q 2 ; t = 0). 

7. Inclusive charged particle cross section for particle with |r^ a {,| < 1.5. The 
model predictions are shown as full line and compared to HI data (see text). 

8. Inclusive transverse energy distribution of jets calculated for jets with \r]j et \ < 
1.5 and compared to HI data. 

9. Inclusive pseudorapidity distribution of jets in rapidity gap events with 
E t j et > 4 GeV calculated with the model and shown with HI data. 

10. Ratio of hard diffractive cross sections for two-jet events and total diffrac- 



tive cross sections. According to forthcoming the UA8 data |47| , f46| to 
whom our calculation might be compared, a lower cut in transverse energy 
of 8 GeV has been applied. 

11. Distribution of the longitudinal momenta of jets normalized to their max- 
imum value in the CM system of the Pomeron-proton interaction. The 
predictions of the model (labeled "MC") are shown for a momentum frac- 
tion of the quasi-elastically scattered proton between 0.92 and 0.94. The 
uncorrected UA8-data || are plotted separately (labeled "UA8"). 

12. Pseudorapidity distributions of jets in the CM system of the proton-antiproton 
interaction. The Monte Carlo results (MC) are shown for a momentum frac- 
tion of the quasi-elastically scattered proton between 0.92 and 0.94. The 
uncorrected UA8-data [[| are plotted separately. 

13. Distributions of the total longitudinal momentum of a two-jet system nor- 
malized to the maximum value for the quark and gluon contribution of the 
Pomeron and both contributions together are shown. Separately we plot 
the uncorrected UA8-data [El. 
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